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Abstract 

The radiative capture process np d'j is considered within the framework of a recently developed 
six-quark dressed-bag model for the nucleon-nucleon interaction. The calculations presented here 
include both the nucleon current and the meson-exchange current contributions. The latter uses 
short-range hadronic form factors for the pion exchange currents consistent with the soft cut-off 
parameter AttAtat from the A^A^-potential. Contributions of the pion exchange current and A-isobar 
current to the total cross section still cannot explain the discrepancy between the theoretical and 
experimental cross sections. Possibilities for new types of meson exchange currents associated with 
chiral fields inside multi-quark dressed-bag states in nuclei are discussed. 
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INTRODUCTION AND MOTIVATION 



Radiative capture of thermal neutrons by hydrogen, ^H(?t,, 7)^H, is now a well-studied re- 
action which is of prime importance in both astrophysical applications and in fundamental 
nuclear physics. This process is one of the simplest classical nuclear interactions involving 
an electromagnetic (e.-m.) "probe", and thus it is doubly interesting due to its close rela- 
tionship to 180° inelastic electron-deuteron scattering near threshold, as well as to deuteron 
photodisintegration 70? — >■ np. The primary interaction is of e.-m. nature, and thus it is 
well understood; the leading effects of strong interaction occurring between two nucleons are 
also well characterised. Therefore, observables for this reaction, if measured with sufficient 
precision, can provide a sensitive test to phenomena associated with, for example, the role 
of sub-nucleonic degrees of freedom. 

The numerous measurements done in previous years have demonstrated that the accepted 
values for cr„p(^H) are free of systematic errors and without masking due to competing 
processes. The modern accepted value of the np d'j cross section for thermal neutron 
capture Q at a neutron velocity of 2200 m/sec is 



np 



;^H) = 334.2 ±0.5 mb. (1) 



We will compare this experimental result with theory. 

For many years a 10% discrepancy has existed between experiment and theoretical calcu- 
lations for the total radiative thermal np capture cross section, np d'-f. The discrepancy 
was first pointed out by Austern [3] as early as in 1953 and has been repeatedly verified in 
the literature since then. The successful explanation of the cross section in terms of meson 
exchange currents (MEC) given three decades ago by Riska and Brown [3] was considered as 
one of the cornerstones of mesonic degrees of freedoms in nuclear physics. Using a realistic 
hard-core wave function for the deuteron, these authors computed the two diagrams with 
one-pion exchange (OPE) initially suggested in 1947 by Villars ^ plus u and A resonance 
diagrams. Although suspected since the Yukawa force was introduced, the work of Riska 
and Brown was the first evidence for the explicit role of mesons, in particular that of pi- 
ous, in nuclear interactions. In recent years the status of meson-exchange currents in the 
deuteron, including the np d'y reaction, has been discussed exhaustively by Arenhovel , 
Mathiot |6|] and Riska 17|]. 

However, apart from quantitative calculations showing the sensitivity of the cross section 
to the A^A^ wave function (at a level of 1.5%) 0, short-range contributions to exchange 
currents {ttNN form factor, p-exchange currents) have been studied only recently. These 
contributions appear to be of some importance, especially in the light of current pictures of 
the quark substructure of the nucleon. 

In the framework of chiral perturbation theory (xPT) it has been shown that the short- 
range interaction is supp ressed in the exchange current contributions From this the con- 
clusion was reached [lO| , that only soft one-pion exchange terms contribute to the two-body 
current. This means that all short-range hadronic processes and higher order corrections 
are cancelled out to rather high accuracy for the magnetic transitions leaving only the pion 
degree of freedom — "a chiral filter'\ This important conclusion shifted the problem of sen- 
sitivity of all the MEC-contributions to the values of cut-off parameters A^^ttv, A^ArTv, Ajtata 
etc. used in current theories for MEC. In fact, all modern one-boson exchange (OBE)-models 
suggested up to date for the description of the A^A^-interaction employ rather high cut-off 



values AjrATAT — 1.5 GeV [ll| and thus the currently accepted approaches for MEC also 
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include similarly high values of cut-off parameters. Meanwhile, all consistent theories for 
n — N interactions and also numerous experimental data result unambiguously in soft cut-off 
parameters around 600 MeV/c However, these soft cut-off values cannot be matched 
with the quantitative description of the iVA^-interaction within the OBE- model [H, [13. 



Moreover these soft cut-offs, being in agreement with the above chiral filter ideas, lead un- 
avoidably to a significant reduction of all the MEC-contributions to the np d'y process, 
and thus there appears again some visible disagreement between the experimental data and 
current theoretical estimations for MEC-contributions. 

In the present work we try to solve this puzzle by use of a new A^A^ microscopic force 
model developed recently jointly by the Moscow-Tiibingen (MT) group 0,^^. Contrary 
to meson-exchange models, the new model includes soft cut-off parameter values only, and 
thus it is in agreement with the above "chiral filter" ideas while providing a quantitative 
and accurate description of iViV-data [3] • 

Another specific feature of the new force model is the fact that mesonic and quark degrees 
of freedom are included into the model in a consistent and explicit form, so that the total 
wave function of the NN system is described as a Fock row: 



V : / 



(2) 



here the second and higher components are of quark-mesonic nature and are essential at short 
and intermediate (r < 1 fm) ranges. On the other hand, it would be very interesting to 
check whether there is a non-negligible contribution to the np —>■ dj reaction at low energies 
from the region r < 1 fm, especially for meson exchange- current contributions. Another 
interesting problem to study here is the role of "inner" tensor force introduced within the 
new NN-mo del [illll. This specific short-range tensor force is based on the symmetry of 



the six-quark wave functions and leads also to some enhancement of the D-wave function of 
the deuteron in the asymptotic region. This modification should affect the radiative capture 
cross section. 

To answer the above key questions, we present in this paper a consistent ab initio calcu- 
lation for the cross section of the np — >■ ^7 reaction with thermal neutrons using the new 
NN -force model. 

The structure of the article is the following: section 2 is dedicated to the general features 
of np-radiative capture and its relationship to the underlying A^A^-force model. In section 3 
we present a brief description of the new force model. The deuteron structure emerging 
from the new A^A^-model is discussed in section 4. The central section 5 is devoted to 
the formalism for e.-m. currents within the quark-meson force model while the calculational 
results are presented in section 6. In the concluding section 7 we discuss the results obtained 
and outline general perspectives of our new approach. 



GENERAL FEATURES OF np CAPTURE 

The main features of low energy np capture are generally supposed to be well under- 
stood ■ We list here the basic points specific to this process. 
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At threshold the neutron is captured predominantly from an / = state. This is to be 
expected due to dominant S'-wave component of the deuteron wave function. This statement 
was verified by explicit calculations reported by Adler [l^ and Jankus Moreover, it 



was found that the dominant transition proceeds from the np ^Sq initial state, while the ^Si 
state is assumed to be of little importance. Since the deuteron ground state is dominantly 
^5*1, the transition involves a change in total angular momentum of A J = 1 and thus it is a 
Ml transition. 

The transition from the ^ Sn to the ^5*1 component of the deuteron wave function dominates 
the impulse approximation [20|. Since the continuum ^5*0 state is of isovector while the 
deuteron is of isoscalar nature, the dominant ^Sq ^ d transition has an isovector character 
and thus the transition involves an isospin-flip amplitude; the emitted photon has / = 1. The 
isovector magnetic moment for the two-nucleon system, (/ip — fin), can be used to simplify 
the calculations. Similarly the ^5*1 — > d transition is isoscalar and involves the isoscalar 
magnetic moment, (yUp + and isoscalar charge. 

The problem under question is closely related to both deuteron photo- and electro- 
disintegration. Indeed, the magnetic photo-disintegration is merely the time-reversed process 
to the np — > d'y reaction and thus involves precisely the same matrix elements and transi- 
tion currents as the latter reaction. The deuteron electro-disintegration requires a broader 
knowledge of the transition current process; more precisely, the relevant magnetic transition 
current involves a value of momentum-transfer squared, g'^g^ = g^. The radiative capture 
process is accompanied by emission of a 7-quant with = c<j^ — = where |q| equals the 
deuteron binding energy, |q| = e^. Thus, the current we are dealing with here is merely a 
special case of the currents involved in electro-disintegration. 

The expression for the cross section can be divided into the phase space and kinematical 
factors and into the current transition matrix element which describes the dynamics: 

da -^AEl^ 1 1 MpM„M, |2,3„ . 

where the transition amplitude takes the form: 

-^/i = (^rf|/(q) -eAl^pn) (4) 

with j(q) for the hadronic e.-m. current and stands for the photon polarisation vector. 
ipd and ippn are the initial and final state wave functions. 



THE EFFECTIVE A^iV-INTERACTION IN THE DRESSED-BAG STATE MODEL 

The long-range part of the A^A^-potential is certainly the most accurately known and it 
is generated from the exchange of vr-mesons. This part of the potential provides a rather 
good description of the static properties of the deuteron, and gives rise to a noticeable D- 
wave component in the deuteron wave function. The structure of the A^A^-interaction at 
intermediate and especially at short distances is still not well understood. In conventional 
models it is assumed to be generated from two-and three-pion exchange, although such a 
construction is very complicated. Moreover, this part is fundamental in the treatment of 
the nucleon structure since the vr-mesons are regarded as the Goldstone bosons associated 
with spontaneous breaking of chiral symmetry. A complete account of the physics related 
to the pion degrees of freedom can be found in ref. 
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The most natural way to take into account the multiple vr-exchange is to use disper- 
sion relation techniques. The important ingredients are the correlated and uncorrelated 
27r-exchange contributions. The part which comes from uncorrelated 27r-exchange involves 
isobar components in the intermediate states; both crossed and uncrossed diagrams con- 
tribute. Another part comes from pion rescattering and is included through the vrvr scatter- 
ing amplitude to the t-channel. It was assumed in previous years that the 27r-exchange 
amplitude with vr— vr resonance-like interaction in the s-channel and A-isobar intermediate 
states can be parametrised, in first order, by the exchange of heavier scalar-isoscalar mesons. 
This gave rise to the concept of the scalar-isoscalar a meson exchange. The A^A^-potential 
constructed in this way is momentum- (and/or energy-) dependent. This is certainly not 
very surprising since at very short distances the potential has no reason to be local. In 
the literature we can find several modern parametrisations of the A^A^-potential (see e.g. 
refs. inililHIi^); they may be classified according to their description of the medium 



and short-range part. However, it has been demonstrated independently by a few groups 
that the exchange of a correlated 27r pair in the scalar-isoscalar channel cannot give any 
significant attraction between two nucleons sufficient to couple them to a deuteron or to 
produce nuclear binding (see e.g. ref. ^J). This failure has perhaps a fundamental charac- 
ter for our understanding of the nuclear force, and thus could mean some deep revision of 
the force at intermediate distances. 

Another fundamental problem with the current status of the force is the short-range 
behaviour of the A^A^-interaction in terms of meson- nucleon form factors. In fact, at short 
distances (r < 0.8 fm), all these potentials are regularised in either a completely phenomeno- 
logical way (e.g. as in the case of the Paris A^A^-potential) or by phenomenological form 
factors at each meson- nucleon vertex with cut-off parameters fitted to A^A^-phase shifts . 
This short-range part of the A^A^-potential is indeed very crucial in nuclear physics since 
it gives rise to a strong repulsion between the nucleons and leads in particular to the sat- 
uration of nuclear matter in the usual non-relativistic many-body formalism. However the 
cut-off parameter values A^tva^, ^pNN, ^ttNA etc., accepted in current OBE-models for A^A^- 
potentials, e.g. A^^vAf — 1.5 — 1.7 GeV, exceed the respective cut-off parameter values derived 
independently from all microscopic and dynamical models for nN interaction, as well as the 
values inferred from experiment, A™^^ ^ 0.5 — 0.7 GeV. 

There is another serious problem with the short-range A^A^-interaction which is attributed 
mainly to a;- meson exchange. The ujNN coupling constant fi'^TVAf/^^; providing a good fit 
to A^A^-scattering data, i.e. qI^mnI^'^ = 13 -i- 15, exceeds three times the value dictated by 
5'[/(3)-symmetry jl^], while all remaining couplings are still in agreement with this symme- 
try. Moreover, Feshbach has demonstrated that the SU{?)) symmetry should be a rather 
good one in this region (see the detailed discussion of the above mentioned disagreements 
and inconsistencies in OBE-models in ref. These strong disagreements peculiar to the 

current OBE-models require a search for another model for the short-range behaviour of A^A^- 
potentials. Such an alternative model for the A^A^-force has been developed recently 3, ^ 
by the Moscow-Tiibingen group. 

In conventional models the A^A^-potential appears when the mesonic degrees of freedom 
are eliminated explicitly from the nuclear wave function. Contrary to this, in the new model 
the main part of the interaction appears when mesonic degrees of freedom are included 
explicitly. The key element of the new approach is a dense six-quark droplet surrounded by 
strong meson-fields. Fig. ^ Naturally, such a bag has been called as Dressed Bag while the 
respective driving mechanism for the NN-foice in the model is the Dressed-Bag-State (DBS) 
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FIG. 1: The graph illustrates two sequential Tr-meson emission and absorptions via an intermediate ct- meson 
and the generation of a six-quark bag. 

mechanism jl5l.[l6l|. Contrary to the conventional OBE-models, here the mesonic degrees of 
freedom are included in an explicit form. Thus the total wave function of the NN system 
should be presented in form of eq. (j21). The A^A^ wave function cam be generalised further 
with respect to the first version of the model of ref. by the important configurations 
which incorporate charged 27r-clouds around the six-quark core (see Figs. and[2|D). 

The channel couplings between different components are calculated from a quark-meson 
microscopic model [16]. Typical diagrams representing the above couplings have a structure 
shown in Figs.Et andEb, where the a- field represents the scalar-isoscalar channel, while the 
s-channel 27r-exchange (shown in Fig. Eb) is dominated by the isovector combination of two 
pions. 

The t-channel (external) tt- and 27r-exchanges are assumed to be usual Yukawa-type 
interactions with soft Attatat cut-offs, which represent the peripheral parts of the A^A^- force. 
These soft cut-off form factors provide very natural separation between t-channel Yukawa- 
like (i.e. peripheral) vr- and 27r-exchanges and s-channel a-, p- and 27r-exchanges localised 
in the bag region, i.e. at r < 1 fm. The attractive interaction in the A^A^ channel in the 
DBS approach comes mainly from the a-field (or cr-exchange in the s-channel) generated 
in the intermediate dense six-quark bag. The dominating dressing in even A^A^-partial 
waves {S, D,- ■ ■) is the cr-meson field while in odd partial waves this will come from the 
p-meson field. It is also important that the initial and final A^A^-components (in S- and D- 
waves) are described by mixed symmetry six-quark states |sy[42]a;;L = 0,2, ST) whereas 
the intermediate bag has a completely symmetric structure |s^[6]a;;L = 0,ST); both six- 
quark components being orthogonal to each other. Thus, the effective interaction in the 
A^A^ channel will result after elimination of other Fock-components in the row (eq. ^ except 
"^NN- The effective one-channel Hamiltonian for the A^A^-component takes the form jlslEq: 



OPE 



ort'- 



(5) 



Here Vqpe and Vtpe are the tt- and 27r-exchange potentials in the t-channel taken with 
the soft cut-off parameter Attatat ~ 0.5 GeV. VNqN is the main part of the interaction 
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FIG. 2: Two diagrams for A^A^ coupling to a six-quark bag with one a- or p-meson (two correlated 
pions) or with two uncorrelated pions. 



at intermediate ranges, which is generated due to the intermediate DBS-production with 
dominating cr-meson dressing of the bag and includes inevitably an energy dependence which 
is the consequence of the elimination of some channels from the many-channel problem. The 
orthogonality pseudo-potential Xort^ with Xoh — > oo (in practice the value of A is usually 
taken about 10*^ MeV- fm'^) and 

r = |0s)(0s| (6) 

is the projection operator onto the |Os)-state in the A^A^ channel, where the F operator 
excludes the |s^[6]a;) six-quark component from the initial A^A^ channel. As form factors |0s) 
we use the h.o. function with the common radius tq: 



|0s) = 



(7) 



The A^A^-interaction in the fully symmetric six-quark state [s^]x (which is undressed here) 
is strongly repulsive [24] both in One-Goldstone-Boson (OGB)-exchange (isl . 26l | and in One- 
Gluon- Exchange (OGE)-models ji?, H, (for q — q interaction), and thus we replace this 
highly non-local A^A^-repulsion by the repulsive pseudo-potential Xort^ in a separable form. 
The forms for Vqpe and Vtpe adopted in this work are given in the Appendix to this paper. 

The effective A^A^-interaction VNqN induced by the intermediate DBS generation has the 
following form: 

VN,N = XlUms){M. (8) 
for the singlet channel. For the triplet ^Si — ^Di coupled channels we have 



X\,,{E)\2s){2s\ AU(i?)|2s)(2rf| 



V 



NqN 



X\^^{E)\2d){2s\ X\^^{E)\2d){2d\ 



\ . 



(9) 



The bag radius is taken as 6 = 0.5 fm and relates to h.o. radius tq = y2/3 6. The standard 
normalised h.o. wave functions in the A^A^ relative motion |2s) and \2d) are formulated as 
follows: 



|2s) = 
\2d) 



ATTilTrn 




(10) 

(11) 
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FIG. 3: NN-phase shifts in the DBS model up to 1 GeV 

The energy dependent coupling constants A^^^,(i?) in eqs. (jH)) and © are calculated from 
microscopic 6g- and ^Po-Quark pion coupling models jl^ and then are compared with those 
found from a fit to A^A^ phase shifts, Fig. |2l from zero energy up to 1 GeV. An impressive 
agreement for all relative values XsLi'^i^) was found. The energy dependence for \gi^j^i{E) 
has been parameterised in ref. in the Pade-approximant form 

>^-kAE) = ■ (12) 

However, the energy dependence \g]^i^t{E) derived from the fit occurred to be somewhat 
stronger compared to the theoretical predictions. In the present work we use the values and 
the energy dependencies for \gj;^^i[E) as tuned from the fit. 

The parameters for -ffe// in eq. (0), which provide a very good fit to A^A^-phase shifts 
from zero energy until as high as 1 GeV for ^5*1 — ^ Di and ^^o-channels, are given in Table HI 

It should be stressed that the range and strength parameter values displayed in Table H] 
are not only in good agreement with the values derived from the quark-meson microscopical 
calculations and empirical high-energy A^A^ phase shifts, but the same values of the param- 
eters result in a perfect agreement for effective range parameters both in triplet and singlet 
channels; they also yield very good predictions for the deuteron structure (see below). 
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TABLE I: The potential parameters found by fitting the A^A*" scattering phase shifts up to 1 GeV 
(Ref. [3). 





ro 
fm 


r2 
fm 


Aoo 
MeV 


1,1 

'^122 

MeV 


^102 

MeV 


Eo 
MeV 


a 


MeV 


/3 
fm-2 




0.41356 


0.59423 


-328.55 


-15.65 


-44.06 


693 


-0.05 


-4.0573 


0.5301 


'So 


0.430 




-328.9 






693 


-0.05 


-8.803 


0.6441 



THE DEUTERON STRUCTURE IN THE DBS MODEL 

For this work on the deuteron we restrict ourselves to two components of the full Fock 
column eq. (0) , which can be represented by: 

The additional configurations generate new currents, which will be presented in another 
paper, however they are of no relevance here. According to eq. ()13|) . the deuteron can be 
found in two different phases: 

1) in the A^A^ cluster-like phase, in which its state is described in terms of the A^A^- 
relative motion variable r, total spin Snn = S = 1, isospin Tjy^ = T = 0, and the 
total angular momentum J^n = J = 

2) in the 6g + cr quark-meson phase, where one must use the dressed bag variables, i.e. 
the quark coordinates ri, r2, ... rg together with their spins, colours, isospins, the 
total spin of the six-quark system S^g, and also the momentum of the cr-meson k and 
the total angular momentum of the bag Jb = J = 1. 

The first component describes the arrangement of the six quarks into nucleon clusters with 
the wave function given by the resonating group method (RGM) js^: 

■^NN = A [V^^(123)^^(456)xiviv(r)] (14) 
where r = (ri + r2 + r3)/3 — (r4 + + rQ)/3. The ip^{i,j, k) is the nucleon wave function: 
^(1,2,3) = <^(fi,6)| [l']cSsg,m]cs)T3g:[l']CST) (15) 
and the relative motion function Xnn{t^) has the following angular decomposition: 

u{r) M n.^ w{r) Q , . 

XNN{r) = ^3^1100 (r) + ^—3^1120 (r) (16) 
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with 



Ml, Ms 



(17) 



and IS* Ms) and \TMt) are the spin and isospin functions for the two-nucleon deuteron state. 

The DBS-component in simplest approximation (i.e. no additional pion loops are con- 
sidered) can be described by the Yukawa-like wave function 

= ^^6g ■ x{ra) (18) 

where M is the normalising factor. We use the abbreviations \(pQq) = \s^[6]x; L = 0,ST) 

and x(^<t) ~ — f^(rTa) ^ ^j^^Yi being the distance of the cr-meson from the six-quark 

core. The a mass is not fixed strictly in this model because it represents a quasi-particle 
defining the quantum numbers of the scalar-isoscalar channel. It depends slightly upon the 
particular system and the type of the process under investigation. Thus, yUo- ~ 400 MeV IsT 
is a sensible choice. 



TABLE II: Deuteron properties in different models. 



Model 


Ed (MeV) 


Pd{%) 


r-m (fm) 


Qd(fm2) 


Md(MJv) 


As(fm-i/2) 


D/S 


RSC 


2.22461 


6.47 


1.957 


0.280 


0.8429 


0.8776 


0.0262 


AV18 


2.2245 


5.76 


1.967 


0.270 


0.8521 


0.8850 


0.0256 


Moscow 99"' 


2.22452 


5.52 


1.966 


0.272 


0.8483 


0.8844 


0.0255 


Bonn 2001 


2.22458 


4.85 


1.966 


0.270 


0.8521 


0.8846 


0.0256 


DBS model 
(A'' A'' component only) 


2.22454 


5.42 


2.004 


0.286 


0.8489 


0.9031 


0.0259 


DBS model 
NN + 6g(3.66%) 


2.22454 


5.22 


1.972 


0.275 


0.8548 


0.8864 


0.0264 


Experiment 


2.2245475(9)'') 




1.971(3)'=' 


0.2859(3)'*' 0.857406(1)"=' 


0.8846(44)^' 


0.0264S' 



") Ref. 0; Ref. |3^; Ref. 0; ^Refs. 0,0; ^) Ref. 0; f) Ref. 0; 

s) Ref. 



Having chosen all the parameters of the effective Hamiltonian Tie// in the A^A^ channel, 
the deuteron two-phase structure can be derived with no additional parameters. In Ta- 
ble inithe main static properties of the deuteron as emerging from the above A^A^-model are 
compared to experiment and to predictions of other modern NN-force models. 

The weight of the DBS-component in the deuteron wave function, as calculated from 
the above model, is 3.66%^]. We emphasise here that the six-quark component weight 
in this model is a direct consequence of the six-quark model parameters established either 
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FIG. 4: Deuteron wave functions u{r) and w{r) calculated in the six-quark dressed bag model. 

theoretically or from the fit to NN phase shifts. And thus this weight cannot be varied or 
fitted to other data. This weight looks to be quite significant and it is within the limits 
found with other hybrid models. A reasonable variation of the radius parameter tq for 
the 6g-core in connection with the readjustment of the strength A does not result in any 
significant change of the weight for the DBS component in the deuteron. This holds also for 
other observables. The level of agreement between theoretical predictions of this model and 
the most accurate experimental data for deuteron observables is not worse than for the best 
modern A^A^-potentials of conventional type. Our results for the np — > d'y process presented 
in the paper can be fully compared with those of other modern A^A^-models. 

It is also of interest to compare the form of the deuteron wave functions predicted by 
our model and traditional ones; see Fig. E] for both the S- and the D-wave components. It 
is evident from Fig. H] that the DBS model generates a deuteron S-wave function with a 
node at short distances (r ~ 0.5 -h 0.6 fm). This node is a consequence of the requirement of 
orthogonality between all S- (and P-) wave functions in the NN channel and the quark part 
of the DBS-staXe, projected onto the NN channel. There is also a rather small short-range 
loop in the deuteron Z^-wave component of DBS-wave function. 

Thus, two-component wave functions are present in our approach, both in the initial and 
in the final (i.e. deuteron) channel. It should be stressed, that both components have a 
different structure and both will contribute to the np — *■ dj transition amplitude, and the 
contribution of the DBS-component will by no means be negligible. Our main concern in 
this paper is the comparison between predictions of our model and the conventional A^A^- 
potential models, with the inclusion of both the nucleon currents and MEC-contributions. 
Because the conventional models have no other contributions we postpone the study of the 
new currents inside the DBS-model to future work. We present in this paper the comparison 
with traditional OBE-models in the A^A^-sector and the conventional MEC-operators. 

The wave functions of the np continuum and the deuteron are obtained by solving the 
Schrodinger equation with Tie// from eq. ©. In the present calculation the same potential 
model is used to determine both the initial scattering state as well as the deuteron wave 
function. 
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GENERAL TRANSITION ELECTROMAGNETIC CURRENTS 



From a quantitative point of view, the relationship between the nuclear wave function, 
the A^A^-potential and e.-m. operators is expressed via a continuity equation, which ensures 
that the total current is conserved (global gauge invariance) . The e.-m. current operator j 
of a nuclear many-body system is related to the Hamiltonian of the system by the continuity 
equation 

V-i + i[H,p]=0, (19) 

where p is the charge operator of the system. So that once the total system Hamiltonian is 
known, eq. ()19|) gives the current operator j for the meson-exchange-current (MEC) calcu- 
lations. 

The contribution of MEC to e.-m. processes on nuclei constitutes an important mani- 
festation of meson degrees of freedom for charged mesons mediating the strong interactions 
between nucleons in nuclei. They participate in the interaction of an external e.-m. field 
with a nucleus and result formally in the appearance of nuclear currents in the form of 
two- or many-body current operators. So they are intimately linked to the underlying A^A^- 
interaction. However, for a given A^A^-potential there does not exists an a-priori way of 
constructing the appropriate MEC-operator, unless the A^A^-potential is derived from an 
underlying more fundamental field-theoretical framework with explicit incorporation of sub- 
nuclear degrees of freedom. Although the existence of exchange currents associated with 
an A^iV-potential has been acknowledged for a long time, the early realistic potentials, be- 
ing phenomenological to a large extent, prevented thus the consistent consideration of such 
exchange currents [ST^ . 

The situation is complicated by the fact that the e.-m. current operator contains ir- 
reducible two-body exchange current components, which through the continuity equation 
depend on the potential model. Therefore, consistent gauge invariant calculations require 
in general, that the same potential model must be used to generate both the wave functions 
and the exchange current operators. 

The full hadronic current operator has the following matrix form in the DBS model: 

where the diagonal currents 

jr(^) = yA^) + jr^'i^) + yfi^) (21) 

J^'i^) = (22) 

correspond to the interaction of an external e.-m. field with the NN and the Dressed Bag, 
respectively. The third term in eq. (j2ip corresponds to the isobar current (IC) associated 
with the intermediate A(1232) excitation. The transition (non-diagonal) currents j^^^^"^"^ 
or j^^^^^^ describe the interphase transitions induced by the external e.-m. field. 

We consider here only that part of the hadronic e.-m. current operator j^^ given by 
eq. (1211) . which is associated with cluster-like nucleon-nucleon configurations. 
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FIG. 5: One-body e.-m. current in lA which is constructed in the DBS model on the basis of the 
diagrams presented in Fig. [5J Region A corresponds to the generation of the initial A^A^ scattering 
wave function: in the region B, the emission of a 7-quantum by a single nucleon takes place and 
leads to the Ml transition; region C generates the final deuteron wave function. 



Impulse approximation 

The simplest description of nuclei is based on a nonrelativistic many-body theory of inter- 
acting nucleons. Within this framework, the nuclear e.-m. operators are expressed in terms 
of those associated with the individual protons and neutrons — the so-called impulse ap- 
proximation (lA). In this approximation the one-body nuclear e.-m. p{q) and j{q) operators 
are obtained from the covariant single-nucleon current 



u{p 



2M, 



N 



(23) 



where p and ff are the initial and final momenta, respectively, of a nucleon of mass rriN, 
and -Fi(Q^) and -^2(2^) are its Dirac and Pauli form factors taken as a function of the 
four-momentum transfer: 



(24) 



The current is expanded in powers of l/rriN and, including terms up to an order of I/itln, 
the nuclear current operator can be written as a sum of one-body current operators which 
takes the form |4ol |: 



AMn 
i 



AM, 



N 



(25) 



where the first term is the so-called convection current and the second is the magnetisation 
current. 

The Sachs form factors Qe and Qm are related to the Dirac and Pauli form factors in 
eq. (j2ni) via 39]: 

Ge{Q') = F,{Q')-^F2{Q'), (26) 
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(27) 



and are normalised so that 

_ 0) = g^{Q' = 0) = 1, gi,{Q^ = O) = /^^ + /x„ = O.SSO i^n, (28) 

SLiQ'' = 0) = - /x„ = 4.706 (29) 

where /ip and //„ are the magnetic moments of the proton and neutron in terms of the 
nuclear magneton fi^. The superscripts S and V of the Sachs form factors Qe and Qm 
denote isoscalar and isovector combinations of the proton and neutron electric and magnetic 
form factors, respectively. Then the transition matrix elements take the form 



where 



and 



1 u{r) 



\lMs)s\00h 



w{r) 



(30) 
(31) 

(32) 



MlMs 



are ^5*1 and ^Di-deuteron wave function components, respectively. 
For the initial np ^Sq continuum state we use a similar form 



1 z(r) 



47r r 



|00)5|10), 



(33) 



The direct calculation for the impulse approximation (lA) results in the following expres- 
sions for the amplitudes M. 



M.so^ss, = G{i_^ss, \/- Ei^Mjlk lA)n,(q) 



(34) 



/Stt 

V k 

where the functions G{g^_^3Si ^^^So-*^Di taken as follows: 



V 



2Mn 



f u{r)z{r)jo{^)dr 
Jo 2 



(35) 



(36) 



Gi^o^^n, = OUQ' = 0) 



2Mjv 



w(r)^(r)j2(y)c?r 



and with 



(37) 
(38) 
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Mia = Gli^, ^ Y^i^Mjlk lA)Fi.(q) 



(39) 



The lA cross section is thus proportional to the quantity Gf^ 



So^d 



as one can neglect 



the states other than ^5*0 in the entrance channel because of its low energy: 



0) 



2Mn Jo 



M(r)z(r)jo(y) - -^w(r)2;(r)j2(y ; 



dr 



(40) 



The result of the destructive interference between G^a^ and Gj^'^ is a sharp 
minimum in the cross section for the e + d^e + n + p reaction which has not been seen 
experimentally. Notice that at high momenta the transition to the Z)-state of the deuteron 
is as important as the transition to the S'-state. The total matrix element changes its sign 
at about ^ 12/m"^. 

Although the lA has been utilised frequently in investigations of weak and e.-m. interac- 
tions in nuclei, it does not preserve all the salient features characterising the "elementary" 
interaction with free nucleons. For instance, even though the e.-m. current for on-shell 
nucleons is conserved as implied by gauge invariance, the nuclear or hadronic e.-m. currents 
generated by the conventional lA are not conserved in general and are inconsistent with 
gauge invariance 
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To restore current conservation in conventional lA, we should incor- 
porate also MEC and take seriously the argument that only the sum of all the contributions 
shall be conserved. In other words, we have to take seriously the two basic ingredients of 
the lA, viz.: (1) The interaction with the whole nucleus can be approximated by a simple 
sum of the "elementary" interactions with the constituent nucleons, and (2) each elementary 
interaction with the constituent nucleon takes place instantaneously. By "instantaneously" 
it is understood, that the elementary interaction takes place in region B of Fig. El and thus 
is not affected by the presence of the initial and final state interactions. Even if the off-shell 
effects arising from the binding of constituent nucleons by the nuclear potential can be ne- 
glected, the masses for the initial and final states of the nucleon cannot always be treated as 
the same because of the instantaneous absorption of the four-momentum transfer. Current 
conservation is already violated in this elementary interaction if, in spite of such difference 
in the effective masses, we adopt as our basic input the form factors for free nucleons. To 
guarantee that the resultant impulse approximation is in accord with current conservation, 
we have to attempt to restore current conservation for each elementary interaction by gen- 
eralising the on-shell form factors to include the off-shell effects due to, in particular, the 
instantaneous absorption of the four-momentum transfer l4z 



Pion-exchange currents with the soft cut-off parameters 

It is important to note here that both the pion-exchange interaction and the model 
independent pion-exchange current operator may be derived as nonrelativistic limits of the 
corresponding operators obtained from a relativistic chiral symmetric Lagrangian model for 
the meson-nucleon system. The two-body currents generated from the meson degrees of 
freedom can be classified according to their time (charge) and space (current) components, 
and also to their isospin structure (isoscalar or isovector). These contributions are of different 
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FIG. 6: The seagull term contribution in the DBS model. Notations for the regions are the same 
as in Fig. [51 for the one-body current. 

nature, the isovector current is of the same 1/M order as the one-body current, while all other 
transitions appear at the relativistic level and should compete with relativistic corrections 
to the impulse approximation. 

The diagrams for MEC of longest range are presented in Figs. IHl and [7| All transition 
operators representing these exchange-current processes have isovector character, and do 
therefore lead to transition from the deuteron state to the ^5*0 state and not to the a- 
state (3^1 + ^Dx). There is, of course, no limit to the number of possible contributing 
exchange- current processes, but only those of pionic range are really significant, since short- 
range heavy-meson exchange processes are suppressed by the strong short-range repulsion in 
the nucleon-nucleon interaction It means that a isoscalar counterpart of the isovector 
currents, due to p7r7 and loi^'y contributions contribute only a little to the isovector transition 
under consideration and thus have been neglected in the present study. Besides that, the 
class of model-dependent currents between different mesons, i.e. p7r7 (u;7r7), are purely 
transverse and are of short range, which is due to the large p- [iS] mass. Therefore, their 
contribution at low momentum transfer should be very small |39]. It should be emphasized 
also that the contribution related to these operators are rather sensitive to the values of 
cut-off parameters. For the soft cut-off parameters used in the DBS model (A^ ~ 0.7 GeV) 
these contributions will be even less important than in the traditional OBEP models where 
A^ ~ 1.2 GeV is usually employed. 



The Seagull term 



The most important MEC process is due to the seagull current operator J^*^" which is a 
two-body current operator (see Fig. EI): 



jr(A:i,fc2)=i 



NN 



mt 



02(01 -ki) Oi(a2-k2) 



(41) 



ml + kf ml + /cl 

The configuration-space expression [39[ may be obtained from the respective Fourier 
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transform of J^^"'{ki, k2 



which results in: 



NN 



. + -)- {B2{B^ ■ iy-''' + ai(a^2 ■ iY"^-''] (43) 



Ti X T2) [m. 



' ' ' r Arrr 

Thus the transition matrix element is of the form: 

It can be naturally decomposed in two terms: 



(42) 



(44) 



^sea 



So^d 



where the functions Gf 



X M ^ J2i^MjU\lX)Yu{ci) 



So^'^^Si ^^'^ ^i5o^3Di defined as: 



7. 



wNn\ 4 



An I ml 



u{r)z{r)ji{ — )(1 + m^r) — - — dr 



(45) 



(46) 



Gf 



sea 



NN 



-^y^ u;(r)2;(r)ji(y)(l + m^r)— ^ 



-dr 



(47) 



T/ie Pion-in-flight term 

Among all the conventional diagrams, the pion-in-flight term is the most difficult one to 
calculate, since it does involve two-pion propagators (see Fig. [7j). This diagram represents 
the direct coupling of the photon to the exchanged pion. It corresponds to virtual pion 
photoproduction on a single nucleon by the pion-pole amplitude, followed by absorption on 
the second one. The result of the calculation for this pion-in-flight (pif ) term is: 



(ai ■ ki){a2 ■ h) 
{ml + kf){ml + kl) 

And the transition matrix elements can be written as follows: 
consisting of two components: 



So^d 



Gl 



pif 



_1_ /^P«/ 



^ J2{lMjU\lX)Yuici) 



(48) 



(49) 



(50) 
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FIG. 7: The pion-in- flight term contribution in the DBS model. Notations for the regions are the 
same as in Fig. [S] for the one-body e.-m. current. 



where the functions GTg^^^g^ and Gi%j_,3j-,^ take the forms: 
//J 



p2 

IttNN 

An 



2j_ 

3 ml 



' f_ 2d_\ 

^ dr"^ r dr ) ° 

( 5 d 3\1,,,,,' 

+ 3"^ + --r + ^]t2} u[r)z[r)dr 
\ dr"^ r dr / 



(51) 



where 



and 



'flNN\ 2 q 
An I 3 ml 



/ d^ d\ 
ydr"^ r dr J ° 

/ 2d 6 

\ dr'^ r dr 



t2 1 --j=w{r)z{r)dr 



a 2 



a = Jml + -r/2)_ 



(52) 

(53) 
(54) 



In the long- wave limit (limit for small qr) the Bessel functions of the eq. (jSHj) can 

be reduced according to: 



.0 z- 



(2/ + 1)!!' 

so that the eqs. and result in: 



(2/ + 1)!! = 1 • 3 - 5 - ■■ ■ (2/ + 1), 



( fwNN 


\ 2 q 


\ An 


J 3 ml 


fnNN\ 




An J 


3 ml 



<^ (m^r - 2) 

r 



1 \ e~™' 

71 X + 



M(r)z(r)(ir 
- > w{r)z{r)dr 



(55) 

(56) 
(57) 
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The A{1232) -isobar current term 

The process corresponding to the virtual excitation of a A(1232) followed by the AN — *■ 
A^A^ transition interaction (Fig. |HJ is described by the A-current J^. There are many 
approaches to evaluate this effect. We assume A-dominance, and we shall use here the 
static quark model and the sharp resonance approximation. 




I 

I TT 

I 

I 

N : N 



FIG. 8: The A-isobar current diagram. 



In the static quark model the exchange current due to the A resonance is |43l |: 



At' 



N 



ml + ki ml + /cl 

{a^ X X g((T^ x k2) (c?^ x ki) x q(a^ x ^2) 



ml + /co 



ml + kl 



(58) 



here ki and ^2 are the fractional momentum transfers imparted to the first and second 
nucleons, and Qt^nn is the pseudoscalar vrA^A^ coupling constant, which is related to the 
pseudovector vrA^A^ coupling constant f-^NN as: 



nNN 



and 



f2 

JttNN 



2Mn m^ in 
The transition matrix element for the A-exchange current is: 



0.075 



Mi 



So^d 



consisting of two contributions: 



and 



(59) 



(60) 



(61) 
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FIG. 9: The one-pion-exchange potential (OPEP) modified by the hadronic form factor of 
monopole type at each vertex. 

where Gi^p^a^,^ in the long-wave Hmit has the form: 

-.2 n / <ry~i \ 3 



4 

X — 

25 



/ A / A / 1 3 3 \ e-"^-" ^ 
«;(r)2;(r) 1 + + - dr 



(62) 



Effect of short range correlations 



One important feature of the vr-exchange currents which has not been much emphasised 
in previous studies of the radiative neutron capture is the sensitivity of these exchange- 
current effects to the nNN form factor. From a general point of view, the ttNN form factor 
originates due to the modification of the bare vr-exchange current at short distances 
where the finite size of the nucleon (and the pion) together with the vertex renormalisation 
associated with very short-range multi-pion exchange processes (the coupling of the pion to 
heavier systems such as p or Svr) can be of some importance j45|. The former corrections 
are roughly related to the confinement radius of the quarks in the nucleon. 

It is rather difficult to calculate the form factors exactly, as they involve many non- 
perturbative contributions js^. For practical applications they are usually parametrised 
in momentum space by the monopole form, the mass scale A being a parameter chosen in 
order to be consistent with the experiment. All the consistent estimates given up to date 
lead to a cut-off parameter A^jvat lying around 0.5 -r- 0.8 GeV/c. In the same way as the 
e.-m. form factor in MEG operators has been related to the longitudinal part of the charge 
operator from the continuity equation, the hadronic form factors in these operators should 
be connected to the choice of a particular A^A^-potential. 

One takes 

F.NNiq) = (63) 
A^ + 

at each nNN vertex. The values of the cut-off masses are, of course, to be related to the 
cut-offs used in the construction of the A^A^-potential. 

One important point in deriving the expressions for MEG with the hadronic form factors 
is to insure gauge invariance, or in other words, the current conservation. It has been 
shown a long time ago that current conservation is insured, if one takes a vr-exchange model 
for the A^A^-potential together with exchange currents (seagull and mesonic) including an 
appropriate e.-m. form factor of nucleons and pions but without any meson-nucleon form 
factor. The next step is the verification whether insertion of the hadronic form factor at each 
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FIG. 10: The seagull (a) and the pion-in-flight (b,c,d) diagrams modified by hadronic form factors 
of monopole form. The diagrams (b) and (d) are required in order to ensure gauge invariance 



meson-nucleon vertex satisfies this requirement too. The vr-exchange potential is modified 
by a monopole form factor at each vertex in the following way: 



fl 



NN 



'M - ml 



mi 



[a^ ■ q){aj ■ q) ^ 
(p + ml 



(64) 



Apart from the normalisation factor (A^ 



2\2 



this modified potential can be visuahsed 



by the diagram shown in Fig. |H1 where 11 represents a particle of mass with the same 
quantum numbers as the vr-meson (the short-range vertex renormalisation is represented by 
an exchange of a virtual heavy particle, of mass A^). By minimal substitution on this dia- 
gram, it is then straightforward to construct the corresponding exchange currents, modified 
by hadronic form factors and satisfying gauge invariance. These are depicted in Fig. ITUk 
(modified seagull term) and Fig. llOb.c.d (modified mesonic term). 

The seagull current operator is just modified by adding the corresponding factor F^^^(g) 
in the numerator of the expression given in eq. (j3^: 



f2 
■ JttNN 



mt 



0^2 (o^i • ki) Kl - ml _ ai((72 • Al - ml 
ml + kl Al + kf ml + kl Al + kl 



The leading isovector vr-exchange pion-in-flight current operator given in eq. 
modified: 



(65) 
is thus 



3':^Kk2 



fl 



NN 



mt 



(ri X f^Yik^ - k2){Al - ml){a^ ■ k,){a2 ■ k^) 



X 



AI + kl 



(A2 + kl){ml + kl) [ml + kl){ml + kl) 



{ml + k!){Al + kl 



Al + kl 



(66) 



The main difference from the usual derivation of exchange currents without form factors 
is the appearance of the two new contributions, shown in Fig. ITUk and Fig. ITUb . which arise 
in the mesonic term. After the inclusion of the hadronic form factor, the exchange-current 
operators in eqs. and (j^Tj) are modified in the following way: 



-A„r 



[1 + m^^r) 



'1 + m^r) 



A^ 

:i+A^r)-^ 



m: 



(67) 
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A 2 A^" ^^^^ 

where 

A = ^A2 + ig2(i_^2). (69) 

Thus, eqs. (UHl), (^Tj) - (jSH), (jU7j) and (jUHj) define the exchange currents within the 

conventional model to be used in this work. This model has the attractive feature that it 
satisfies current conservation consistent with the one-pion exchange potential. The effect of 
the short-range hadronic form factor to the A-isobar contribution in the A^A^ wave function 
was investigated by transition potential methods ji^ with vr-exchange contributions. For 
the calculations we employed a potential proposed in ref. jisj ]. 



CROSS SECTIONS FOR THE THERMAL NEUTRON CAPTURE 



The thermal np capture is a fundamental process for our purpose because it involves two 
opposite features, long-distance behaviour and short-range processes, where the latter are not 
completely understood. Clearly, a better determination of the asymptotic behaviour of the 
theoretical A^A^ wave function is needed before a quantitative comparison with experiment 
can be made for the thermal np capture cross section. 

For thermal neutrons we can use a nonrelativistic kinematic regime, and we may even 
ignore the deuteron recoil energy. Then the 7-ray emitted in the capture process will have 
an energy essentially equal to the deuteron binding energy In this kinematic regime 

the differential cross section has the form 



|v„ - Vpl (27r)2 2uj^ 



\Mfi\^ (fpd(fq^ 



(70) 



'PI V'V ^"^7 

The integrals over and |q| in the laboratory system where |vp| = 0, are easy to take |4j 
if we ignore the deuteron recoil and set g = a; = e^, i.e. the deuteron binding energy. So, 
these integrations results in: 



daji 



uo. 



7 



dn^ |v„|2(27r)2 l-^^*' 



(71) 



Since we neglect the deuteron recoil momentum, there is no angular dependence in the 
differential cross section dafi/dfl^ due to phase space. Therefore, all the angular dependence 
is due to the dynamics of the capture process. 
The total transition amplitude takes the form: 



Ml 



So- 



Mi 



lA 



Pif 
So^d 



Mi 



So^d 



(72) 



where 



M^ 



So^d 



^So-^d + '-'I 



ea I ^pif 




— Y.i^MjU\lX)Yu 



(73) 
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To compute the unpolarised cross section from eq. ()71|) one has to sum over the final and 
average over the initial spin directions. Carrying out these summations one gets the ex- 
pression for the differential cross section for thermal np capture to the ^Si — ^Di deuteron 
state: 



da 



dVt^ |v„|4(27r) 



So- 



(74) 



After angular integration and introducing the relative np momentum in the cm. system of 
the incident channel: 



p = -Mn |v„ 



1/2 



one arrives at: 



(7{np d'-y) = a 



2p 



GIA I risea 

^So^d + ^1 



So^d 



(75) 



(76) 



where a = e'^/Air is the fine structure constant. 

The relative change of the lA result due to the exchange currents can be expressed in 
terms of the quantity 



(5(MEC) 



I QPif 

^ 5*0— >d ^sea 



r<iA 

'^^So^d 



The empirical value for 5 is normally deduced from the relations: 

a'-P(np ^d-f) = il + 5Y- a^^{np ^7) 

where 



a- 



2p 



Gi 



lA 

So^d 



(77) 



(78) 



(79) 



and G{g^_^^ is determined by eq. pUj) with q = ui^. We have now at our disposal all formulae 
for the analysis of both the lA- and the MEC-contributions in the DBS model. All the above 
formulas for transition matrix elements are derived for the general case and can be applied 
to related processes as well. In the Jong- wave limit, the equations for the MEC will lead to 
the forms familiar from literature 



RESULTS AND DISCUSSION 

The most important uncertainty in lA comes from the asymptotic behaviour of the wave 
function, because the momentum transfer carried by the photon is almost zero. Therefore, 
it is not surprising that the analysis of the theoretical cross sections within lA jsj has 
demonstrated a very high stability of the lA-results for different iViV-potential models (see 
Table Uni). 

In Table irm the lA results and the normalisation 5 to the experiment are compiled for the 
potentials of de-Tourreil-Sprung jiol], Argonne V14 nil], Paris and Reid Soft-Core lU, 



respectively. The small variations among different iViV-model predictions are due to the 
differences in D-wave admixture (which changes the value of parameter As oi asymptotic 
behaviour of the deuteron wave function) and due to the different short-range repulsion. 
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TABLE III: The values of the one-body contribution for the different potentials to the thermal 
neutron capture reaction as compared with experimental data. 



potential 


Ref. 


lA [mb] 


m 


TS 


[49] 


305 


4.68 


Argonne V14 


ri 1 1 


304.1 


4.83 


Paris 


m 


302 


5.2 


RSC 




300 


5.55 


DBS model 


present 


304.3 


4.79 


Experiment 


m 


334.2±0.5 



The potentials are ordered for increasing repulsion from the TS to the RSC potential. Thus, 
the 10% discrepancy between the experiment and the one-body contribution is significant, 
and cannot be accounted for by any change in the nuclear wave function compatible with 
the fit to the accurate iViV-data. 

It is evident from Table IIIII that the lA-results for all realistic A^A^-potential models 
and the DBS-model as well agree to each other within some narrow corridor. It should be 
expected because the long-range behaviour of A^A^ wave functions for all the models is quite 
similar. However, the close similarity of the lA results for all A^A^-models considered here is 
still informative because the DBS wave functions in both initial and final channels have a 
node, and thus the short-distance behaviour is sharply distinct from those for the traditional 
A^A^-potential models. Thus the 30 mb residual of the cross section should come from MEC- 
contributions for all A^A^-models, no matter whether they are based on meson-exchange 
mechanism, hybrid six-quark model or on the dressed-bag intermediate state production. 
This is our first conclusion. 

The 5- value of eq. (f75j) necessary to match the DBS-calculation to the experimental cross 
section of the np —>■ d'j reaction amounts to 

S = 4.79% (80) 

Our results for the MEC-contributions within the DBS-model which include the seagull, 
meson-in-flight and isobar current terms are presented in Table IIVI For clarity of the 
conclusions we show in this Table the results for both soft (A^r = 697 MeV) and hard 
(Att = 1200 MeV) ttNN and ttA^A form factors. It is evident that going to hard cut-off 
parameters increases all the MEC-contributions quite visibly, as was expected before. The 
MEC and IC contributions in our calculation with the soft A^^^n give in sum 

5(MEC + IC) ^ 2.8% (81) 

and hence there remains some 2% of the total cross section still unexplained with the present 
A^A^-model. Therefore it is appropriate here to discuss which further contributions should 
appear in our approach. 
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TABLE IV: Relative contributions of seagull, mesonic and model dependent resonance terms for 
pion exchange currents in the DBS model for two cut-off factors A,r- The values are given for 
5 — > 5 and S ^ D transitions. The 6'^ present the values for point-like nucleons. 



Term 




A^ = 697 MeV 


5^1^^" (%) 

A^ = 1200 MeV 


Seagull 
Mesonic 

Total MEC 


2.707 0.459 
-1.243 0.459 

1.464 0.918 


2.120 0.393 
-1.065 0.393 

1.005 0.786 


2.422 0.447 
-1.117 0.447 

1.305 0.954 


Isobar current 


1.151 


0.974 


1.062 


Total MEC+IC 


1.464 2.069 


1.005 1.760 


1.305 2.016 



The main difference between the DBS-model predictions and the conventional A^A^- 
potential models comes from the mixed two-component quark-meson structure of DBS- 
approach. The most evident specific contributions are 

i) The first non-vanishing contribution, which is solely due to the two-component struc- 
ture, are transitions between dressed-bag states themselves (these currents can be iden- 
tified with the operator j^^ in eq. (j221))- These can be illustrated for ^5*0 {^Si — ^ Di) 
transitions by the following graphs in Fig. ^2 This process corresponds to the Ml 
single-quark transition in the bag. We would like to stress the point that the spin-flip 
of the quark in the bag is a non-exotic process because in the conventional quark- 
diquark model for the nucleon, the usual nucleon spin-flip (which is responsible for the 
conventional (a ■ V)- coupling in t^N- amplitude) originates also from the unpaired 
quark spin-flip transition. 

ii) The dressed six-quark bag in our model has mainly an uncharged a-meson cloud 
with the effective mass of the meson being m„ ~ 400 MeV. The charged six-quark 
core and the chargeless cx-cloud are oscillating in the intermediate DBS-state around 
the total center-of-mass. Thus these oscillations of the charged six-quark core will 
induce some e.m. radiation similar to radiation emitted by a charged oscillator in 
conventional quantum theory. Therefore this radiation from the DBS-component in 
the A^A^ system may have monopole and also higher multipole character. A Ml 
transition should be especially appropriate to describe polarisation degrees of freedom 
in the n p — >■ c?7 reaction. The point here is, that the scattering and bound state 
wave functions in the A^A^ channel in our model are non-orthogonal to each other due 
to the non-hermicity of the underlying A^A^- "potential" VNqN{E). However, the total 
two- (or many-) component functions '^{NN) should still be orthogonal to each other 
and thus monopole i?0 transitions are quite possible to occur separately in the A^A^- or 
in the DBS-channel while the monopole transition between the initial and flnal state 
vectors should vanish. The above monopole transition would include the transition 
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FIG. 11: The 6-quark contributions to the DBS e.-m. current induced by transitions between 
dressed-bag states. 



^5*1 -^^ Si in the deuteron. However, the isoscalar Ml-transition from the initial ^5*1 
channel to the final ^5*1 channel inside the bag, which is related to the spin rotation of 
the bag as the whole from Ms = +1 (or -1) to Mg = (or vice versa), is still feasible. 

iii) The dressing of the bag includes in our model not only neutral a, but also charged vr- 
and p-meson clouds. It would be quite possible for incoming 7-quanta to interact first 
with these charged clouds of the bag. These specific new meson-exchange currents will 
contribute to the np — > d'y reaction and also to many other e.m. processes 5^, 5^, H^ ]. 



We postpone the discussion of these new effects to our further publications. 

In summary, we presented in this paper the first e.-m. calculations for np — > d'y cross 
sections on the basis of the new A^A^-interaction model. The model employs the soft meson- 
nucleon cut-off parameter values so that we still observe some underestimation of the ex- 
perimental cross section. However, there are a few specific new contributions in the model 
which are fully absent in conventional A^A^-potential models. Our further studies point in 
this direction. 
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ONE- AND TWO-PION EXCHANGE POTENTIALS IN THE DBS-MODEL 



If isospin-breaking terms are ignored and staying within the lowest order of perturbation 
theory, the OPE A^A^- potential in momentum space takes the form: 

V^'i^l = -%i^^^i^(r. . r,), (82) 

where is the mass of the exchanged pion and /^^vAf the coupling constant. 
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When the potential is Fourier transformed to configuration space, it is usually regularised 
to remove the singularities at the origin. This can be achieved by introducing a form factor 
F{(f) which has a dipole form: 



F{q 

where A^attv is the cut-off parameter. 
A typical Fourier transform reads as 



'At 



■NN 



m: 



A2 4- 



(83) 



^tqr 



So, in configuration space the OPE potential takes the form: 



47r 3 



[fc{r)ai ■ aj + fT{r)Sij] Ti ■ t. 



where the tensor interaction operator is 

Sij = 3{ai ■ r){aj ■ r) - {ai ■ aj) 
and the central interaction term is 



fc{r) 



X 



X 



2 - 1 



-a 



with 



X = m^r and a 



A, 



■NN 



(84) 



(85) 



(86) 



(87) 



(88) 



Using the definition ()84|) . the Fourier transform for the tensor potential can be expressed 
simply in terms of derivatives of the central potential, i.e.. 



/T(r) 



with 



fcir) 



1 d fl 
mi dr \rdr 



— e 



fcir) , 



X 



X 



a 



2a 



(89) 



(90) 



and the tensor interaction term can be written in terms of the variables x and a as 



Mr) = — 1 + - + -7 



X x^ 



a 



ax 



1 + ^ 



ax {axy 



2 



1 + — 

ax 



(91) 



We use the averaged pion mass in the OPE potential: 

= (m^o + 2m^±)/3 = 138.0363 MeV, 
and the respective "weighted" value for the pion-nucleon coupling constant 

f^^^/{An) = 0.075 
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The value of the regularisation parameter a amounts to 

a = 5.054. 

The inclusion of the peripheral two-pion exchange contribution Vtpe in our case can be 
imitated by the simplified form [l6l |: 

Vtpe = VSp^iPr^e-^^' (92) 

and gives only a quite small contribution (about 2-3 MeV) in the intermediate region r ~ 
1.5 -h 2 fm which is well beyond the DBS radius. However, this small attractive contribution 
is quite important for the precise description of effective-range parameters and asymptotic 
properties of the deuteron. 
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